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Abstract. The spatial cross-correlation function £c recently introduced by Roukema 
et al. pQ is applied to the equilateral toroidal topology of the universe. Several CMB 
maps based on the WMAP 5yr data are analysed and a small likelihood in favour of a 
torus cell is revealed. The results are compared to ten ACDM simulations which point 
to a high false positive rate of the spatial-correlation-function method such that a firm 
conclusion cannot be drawn. 



o 

m 

6 

£j ■ PACS numbers: 98.80.-k, 98.70.Vc, 98.80.Es 

1. Introduction 

The cosmic microwave background radiation (CMB) provides not only the information 
to constrain the cosmological parameters, but also has the potential to reveal the global 
■ structure of the Universe, i.e. its topology [2, El HI El El Ej- One method to reveal 
the topology is the search for the so-called circles- in-the-sky signature [8]. The CMB 
descends from the surface of last scattering (SLS) which is the spherical surface around 
an observer at which matter is seen to recombine. In the case of a multiply connected 
universe this sphere intersects with its "copies" due to the group of deck transformations 
defining the topology. Since the intersection of two spheres is a circle, one obtains pairs 
of circles along which the temperature fluctuations are identical, or at least, if we could 
receive the CMB as a pure signal proportional to the gravitational field on the SLS. 
Modifications due to the Doppler and the integrated Sachs- Wolfe effect are different 
on two paired circles thus preventing identical temperature fluctuations along paired 
circles. A further complication arises from the uncertainties of foreground emissions 
which additionally influence the observed CMB. Applications of the circles-in-the-sky 
signature to the WMAP data [9] can be found in pH HI M, Q31 Q3] • 

The circles-in-the-sky method utilizes only the one-dimensional pixel information 
along paired circles. To improve the topological signal it is proposed in [1] that instead 
of analyzing correlations along the circles, one should define two spatial correlation 
functions £a and £c which also depend on pixels not lying on paired circles. Then the 
comparison of £a with £c provides the topological signature. The idea is the following: 
The topology is defined by the group of deck transformations T, i. e. a discrete subgroup 
of isometries without fixed points. All spatial points q%,q2 G Ai^ obeying q\ = 7(92) 
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for 7 G T are identified, i.e. the spatial comoving space section M.^ is tesselated by 
considering the quotient Ai^ /V. Define now the set T which contains the elements of T 
but the identity removed. With respect to this T the distance c?topo between two points 
qi , q 2 € M. ^ is defined as 

dtopo^i, <?2) := mmd(qi,j(q 2 )) , (1) 

where d(q, q') is the usual spatial comoving distance in the universal covering space 
fA^. Note, that for sufficiently nearby points q\ and q 2 one has dt opo (qi, 92) > d(qi, q 2 ) 
since the identity is not contained in V. The reverse applies for sufficiently far separated 
points qi and q 2 . Now the two spatial correlation functions £a and £c can be defined as 

U(r) := (5T( qi )6T(q 2 )} with r = d(q u q 2 ) (2) 

for the auto-correlation function and 

£c(r) ■= (5T(q 1 )8T(q 2 )) with r = d top o(<?i, 92) (3) 

for the cross-correlation function. Here, ST(q) denotes the temperature fluctuation of 
the CMB in the direction h corresponding to q and (...) the averaging over the pixels 
with r = d(qi,q 2 ) or r = d topo (qi,q 2 ). In the practical analysis a binning with respect 
to r is necessary due to the pixelized sky maps. 

The auto-correlation function £a(?") is the usual spatial correlation function which 
has a pronounced peak at r = 0. In the cross-correlation function £c( r ) this peak has 
to be produced by points far separated on the sky which are topologically adjacent due 
to G?to P o- One has detected the correct group T when both correlation functions £a(?") 
and £c( r ) are identical, if it were not for the same modifications of the CMB which also 
affect the circles-in-the-sky method. The main advantage over the circles-in-the-sky 
method is that the correlations £a and £c are computed from a much larger set of pixels 
which levels out the modifications to some degree. On the other hand, since the spatial 
quantities £a and £c are three-dimensional measures which are only sampled from the 
two-dimensional SLS, one cannot expect £a and £c to be identical for statistical reasons. 

To illustrate the difference between the two methods Fig. CD displays the distance 

d m m(q) ■= min d topo (q,q f ) (4) 

for the toroidal topology where all three side lengths are either equal to L = 3 or 
L = 4 (in units of the Hubble length Lh = c/H , see section |2]) and the symmetry axis 
pointing to the north pole. Note that the scale in the two panels of Fig.[T] is different. 
In the case of L = 3 the maximal value of d m i n (q) is only 0.45, whereas the case L = 4 
possesses the larger bound 0.72. This demonstrates that the distribution of d min (q) 
extends over an increasingly larger range with increasing values of L since the elements 
7 G T shift a given point over a correspondingly larger distance. The number of pixels 
used to discriminate the topology for a given orientation is much larger for the spatial- 
correlation-function method than for the circles-in-the-sky method which utilizes only 
the pixels along the circles with d m \ n = 0. In this work all pixels inside the KQ75 mask 
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Figure 1. The distance (i m i n (q), eq. (QJ, is shown for the equilateral toroidal universe 
with side lengths L = 3 and L — 4 using the Mollweide projection. The paired circles 
correspond to d m ; n = 0, i. e. the dark blue regions. The black regions show the pixels 
excluded by the KQ75 mask used throughout this paper. 

[HI E] are excluded which are contaminated by the milky way and other foreground 
sources. These excluded domains are shown in black in Fig.[TJ 

The spatial-correlation-function method works best if there is a large number of 
pixel pairs with small topological distances d topo , since the peak of the correlation 
function at r = can then be tested with a good statistical significance. In units of 
the Hubble length one needs the correlation functions for r < 1. For this reason Fig. [2] 
presents the distribution of the number of pixel pairs with respect to d min for three 
different sizes of the toroidal fundamental cell. As noted above, with an increasing 
topological scale L, the maximal value <i m ; n increases since the points are topologically 
connected over ever larger distances. A large number of pixel pairs close to dmm = is 
observed only for sufficiently small tori. Since d m i n is the lower bound of the distance for 
which the data are sampled for the estimation of the correlation functions, this implies 
that the statistical significance for the estimates of £c( r ) for small values of r diminishes 
with increasing L. Thus for too large fundamental cells, the spatial-correlation-function 
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Figure 2. The distribution of the number of pixel pairs in dependence on d m j n is 
shown for three different sizes of the toroidal cell. For the model with side length 
L = 3, 4 and 5 all points on the SLS have a largest distance <i m i n of 0.45, 0.72 and 
1.66, respectively. With decreasing torus size the number of pixel pairs with d m [ n < 0.2 
increases significantly. 

method has the same problems as the circles-in-the-sky method. In the latter case the 
number and the sizes of paired circles declines with increasing side length L until there 
are none at all. For L — 3, 4, and 5 there are 16, 9, and 3 pairs of matched circles. In 
the following the analysis is restricted to L < 5. 

In Roukema et al. p] the spatial-correlation-function method is applied to the 
Poincare dodecahedron and a special orientation is found for which an enhanced 
likelihood occurs. In this paper, the flat toroidal topology with the additional restriction 
that all side lengths are equal, i. e. the cubic topology, is in the focus in order to see how 
this topology behaves using the method of pQ. 

2. The Markov Chain Monte Carlo method for the toroidal universe 

The toroidal topology is defined on the universal covering space A4^ = M 3 on which 
the points q = (x, y, z) G M 3 and q' = (x', y', z') G M. 3 are identified by 

q' = (x + n x L x , y + n y L y , z + n z L z ) with (n x , n y , n z ) G Z 3 . (5) 

Since the group Y is infinite, a finite subset of Y has to be chosen for the computation of 
the distance d t0 po(qi, #2) in eq. ([I]). Here, the subset is restricted to the group elements 
which lead to copies of the fundamental cell directly adjacent over a face, an edge 
or a corner, i.e. to the elements with \n x \ < 1, \n y \ < 1, \n z \ < 1 and excluding 
\ji X i ri yi n z ) = (0,0,0), of course. The minimum in eq. ([T|) is thus determined over 
26 group elements. As already stated, only an equilateral toroidal cell is considered 
in the following, i.e. L = L x = L y = L z . For a general toroidal topology, these 
considerations have to be modified in an obvious way. The lengths L are given in units 
of the Hubble length Lh = c/Hq. For the conversion of the CMB correlation to the 
spatial correlation, the distance to the SLS is required, which in turn is determined by 
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the chosen cosmological parameters. For the ACDM model one gets from Table 2 in [TH] 
the cosmological parameters based on all astronomical observations, see their column "3 
Year + ALL Mean", i. e. fi b = 0.044, fi cdm = 0.223, fi A = 0.733, h = 0.704, n s = 0.947, 
r = 0.073. For these values the distance to the SLS is L SL s = At/L h — 14.2Gpc where 
At/ = t/o — t/sls = 3.329 (rj is the conformal time). For other sets of cosmological 
parameters, one has to rescale the topological length L with respect to the new value 
of An. 

As in [1] the Markov Chain Monte Carlo (MCMC) method with the Metropolis 
algorithm [T7] is used in order to determine the most likely side length L together 
with the orientation of the fundamental cell defined by the three Euler angles 
(a, /3, 7). This leads to a four parameter (L, a, (3, 7) search. The Metropolis algorithm 
requires the specification of a probability estimator P which determines the probability 
min(l, P new /P id) with which a new state is accepted in the Markov chain. In [1] this 
probability is chosen as 



Here, the index i runs over the bins for which the values of £c and £a are sampled and 
Ni denotes the number of pixel pairs contributing to bin i. See ref. [1] for details and a 
discussion motivating this choice. Since <7j is larger for bins having fewer entries N, the 
probability P weights the bins according to their statistical significance. The smaller 
the distance d topo is, the smaller is also the number of pixel pairs contributing to the 
corresponding bin. However, it is just the behaviour at small distances, i. e. the peak at 
r = in ^A( r ), which has to be reproduced by £c( r ) if the supposed group T matches the 
true one. For this reason a second probability P is used in this paper, which emphasizes 
the peak at r = by using instead of the above <jj, eq. (171), the following alternative 



For large values of £a(?0 a small <3"j is obtained thus selecting models producing strong 
correlations at r = 0. This probability is complementary to the former one. It, however, 
ignores the number of entries on which the value of the bin i is based. 

Since the number of entries into a bin increases linearly for an equidistant binning 
in r for an equidistant binning with respect to r 2 is used in this paper. This leads 

to Ni of the same order. For such a binning the complementary behaviour of a and 
a is shown in Fig.0 The sky maps are in the HEALPix format [18] and are analysed 
for Ti s ide = 128. Since these maps have 12n 2 ide pixels, only every 10th pixel is actually 
used in order to keep the numerical effort manageable. Without using a mask this leads 
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Figure 3. The different behaviour of a, eq. ([7|), and tr, eq. (JSj) , is illustrated. The grey 
region shows i a which increases with with decreasing distances r. The reverse 
behaviour is displayed by £a =t & shown as the hatched area. 

to iVpi x = 19661 pixels, and using the KQ75 mask which cuts out 28% of the area, 
-Wpix = 13534 pixels contribute to the computation of £a and £o The computation of 
the spatial correlation function is based on |A^ pix (iVpi x — 1) pixel pairs, which determine 
the increase of the numerical effort with respect to the number of pixels iV P i x . It takes 
1-2 weeks of cpu time on state of the art cpus to generate a Markov chain with 20 000 
states for the above values. In the likelihood ([6]) n = 20 equidistant bins are used from 
r 2 = to r 2 = 0.9. This binning is a compromise in order to resolve in a sufficient 
way the structure of the spatial correlation functions and to get enough entries into 
the bins for a statistical significant estimate. Furthermore, the binning should not be 
finer than the corresponding resolution of the HEALPix map, since distances which 
are not resolved by the sky map cannot be differentiated between adjacent bins. For 
n S ide = 128 the pixel size is 6 P i x = 27.5' [H] which translates into a spatial pixel size of 
Ar ~ 0.0266. Two points within two adjacent pixels can thus have a distance between 
and 2Ar ~ 0.053; thus n ~ 20 represents the maximal useful binning. 

A MCMC run is started by randomly generating an initial point (Lj, fy, 7$) in 
the parameter space. The initial values of Li are uniformly chosen from the interval 
Li G [3.15, 3.45]. By using the Metropolis algorithm, the Markov chain is not restricted 
to this interval, however. The initial Euler angles are uniformly chosen from the 
interval [0, |]. Because of the symmetry of the toroidal fundamental cell, there is 
some redundancy in the sense that the set of Euler angles does not uniquely define the 
orientation, i. e. several sets of Euler angles can give the same values for the Galactic 
coordinates of the symmetry axes of the fundamental cell. In the following the symmetry 
axes denote the directions in which one arrives at the centre of the nearest copies of the 
torus cell by a shift of length L, i. e. the face-to-face neighbours. 
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Figure 4. The distribution of the states of a Markov chain obtained from the ILC 
5yr map with the KQ75 mask is shown in dependence on the length L and the Eulcr 
angle 7. The likelihood function with a, eq. (JT]), has been used. 

3. The spatial-correlation-function signature for the toroidal universe 

Applying the MCMC method described in Section [2] to the cubic topology one observes 
that, after the initial "burn in", most Markov chains occupy a region in the parameter 
space corresponding to a side length L ~ 3.84 using the WMAP-ILC map based on 
the 5yr data [9] outside the KQ75 mask. Fig.H] shows the distribution of a Markov 
chain with respect to the two parameters L and 7 ignoring a and (3. The distribution is 
localized around a small region in the parameter space with Euler angles corresponding 
to the Galactical coordinates of the toroidal symmetry axes given in Table [TJ However, 
this nice behaviour is not observed for all Markov chains. If the initial length Lj is too 
large, i.e. above 4.75, the Markov chains drift towards ever larger values of L. This is 
due to the decreasing statistical significance with increasing size of the fundamental cell. 
As shown in Fig.[2]the number of pixels N r := #{q\d m i n (q) < r}, which give information 
about the small distance behaviour of the cross-correlation function £c( r )> decreases with 
increasing values of L. Thus, the probability that the few remaining pixels generate by 
chance a peak at r = in £c( r ) rises with increasing values of L. This unsatisfactory 
behaviour might be remedied by choosing another probability function as ([6]) which has 
then to depend on L. This topic will be addressed in a future study. 

After having found a preferred region around L = 3.84, ten MCMC chains of length 
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Figure 5. The normalized distribution of the states of ten Markov chains obtained 
from the ILC 5yr map with the KQ75 mask is shown in dependence on the length L. 
The likelihood function with a, eq. ([7]), has been used. The states of ten chains are 
merged so that the shown histogram is based on 200 000 states. These chains lead 
to L = 3.84 ± 0.10. The dashed curve displays the Gaussian parameterized by these 
values. 




Figure 6. The orientations of the symmetry axes of the Markov chains used in Fig. [5] 
are shown on the sky in the Mollweide projection. 
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Table 1. The mean values of the orientation of the torus in Galactic coordinates (I, b) 
obtained from the ten Markov chains shown in Fig. [5] The accuracy is of order 2°. 
The six directions are either opposite or orthogonal to each other. 



20 000 are generated which use initial values (Li, aij, fa, 7j) within the likely domain. 
These ten chains provide together 200 000 states from which the torus orientation and the 
torus size L is estimated. Fig. [5] shows the distribution of L for these Markov chains. An 
estimate of L = 3.84 ± 0.10 is obtained. The directions of the corresponding symmetry 
axes, i. e. the directions to the six centres of the nearest face-to-face neighbours, are 
shown in Fig.[6]in Galactic coordinates using the same Mollweide projection as in Fig.[TJ 
In Fig.[7]these orientations are directly presented in Galactic coordinates which, however, 
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Figure 7. The orientations of the symmetry axes of Fig. [6] are shown in Galactic 
coordinates (l,b). Note that the meridian I = goes vertically through the centre of 
the Mollweide projection in Fig. [6] 
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Figure 8. The distance d m i n (q), eq. (j4|), is shown for the equilateral toroidal universe 
with side length L = 3.84 and the optimal orientation shown in Figs. [5] and [7] using the 
same Mollweide projection as in Fig.[TJ The black regions show the pixels excluded by 
the KQ75 mask. 



distort the maxima towards the poles. Figs.E] and [7] reveal the preferred orientation, 
whose coordinates are given in Table [TJ Fig. [8] displays the function d m i n (q) for this 
orientation as in Fig.[T] masked by KQ75. If one uses instead of the KQ75 mask the less 
restrictive KQ85 mask, which contains more regions with possible residual foregrounds, 
a second slightly shifted orientation arises. But since the application of the KQ75 mask 
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Figure 9. The normalized distribution of the length L is shown according to the 
states of five Markov chains obtained from the foreground reduced maps for each of 
the bands Q, V, and W using the KQ75 mask. The likelihood function with a, eq. J7]), 
has been used. The states of five chains belonging to a given band are merged so that 
the histograms are each based on 100 000 states. These chains consistently lead to 
L Q = 3.83 ± 0.10, L v = 3.85 ± 0.10, and L w = 3.85 ± 0.09 for the Q, V, and W band, 
respectively. 



gives the cleanest pixels [15], the results obtained with this mask are the most reliable. 

In order to test whether the spatial-correlation-function signature is also present in 
the sky maps of a single frequency band, five MCMC runs are applied to each of the 
three bands Q, V, and W using the KQ75 mask. For these three bands, foreground 
reduced maps are available from the WMAP team [9] based on the five year data. It 
turns out that the MCMC chains, having again initial values of Lj randomly drawn from 
the interval Lj 6 [3.15, 3.45], do indeed find the same orientation and the same length L 
for the toroidal topology as the chains obtained from the ILC map. The distributions of 
the chains are localized in a similar way as those obtained from the ILC map for which 
an example is shown in Fig. HI Using initial points in the high probability domain, 
further five MCMC chains with length 20 000 for each of the three bands are generated. 
The distribution of the topological length L derived from these 15 chains is displayed 
in Fig. [U and the obtained values agree well with the result derived from the ILC map. 
This shows that four CMB maps, i. e. the ILC map and the three foreground reduced 
maps for the bands Q, V, and W, give consistent results. 

It is worthwhile to remark that the results do not depend on the chosen probability 
function, i.e. using a or a in (JB]) . Five MCMC chains are generated using a for random 
initial points in the case of the ILC map using the KQ75 mask, and the same scale 
and orientation of the fundamental cell is revealed. Thus the algorithm is robust with 
respect to these two probabilities. Furthermore, the convergence of the Markov chains 
is checked by computing the MCMC power spectrum and the convergence ratio r along 
the lines described in [19] . 

Let us now discuss whether this topology is directly observable for redshifts below 
z = 6. The topological mirror image of our galaxy cluster is beyond the SLS because the 
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Figure 10. The redshift z at which antipodal objects occur in the direction of the 
nearest face-to- face neighbours is shown in dependence on the distance L/2 using the 
cosmological parameters stated in Section [2j 

topological scale L ~ 3.8 is greater than the distance An = 3.329 to the SLS assuming 
the cosmological parameters stated in Section [2j The nearest mirror images of objects 
at the centres of the "faces" of the torus, however, occur at half that length L/2 ~ 1.9 
from our point of view corresponding to a redshift z ~ 5.1, i. e. there should be the same 
objects at the antipodal points of Table [1] at this redshift. The strong dependence of 
this redshift on the distance L/2 is shown in Fig.fTUl 

4. Discussion 

Does the analysis of the preceding Section indeed point to a toroidal topology of the 
Universe or is the enhanced likelihood around L ~ 3.8 a statistical fluke? It could be 
possible that the different contributions to the CMB, i.e. mainly the usual Sachs- Wolfe 
contribution, the Doppler effect, and the integrated Sachs- Wolfe contribution conspire 
in such a way that the cross-correlation function £c( r ) gets by chance relative large 
values at small distances for the obtained orientation and torus size. 

In order to settle this question one has to generate a large number of ACDM 
random sky simulations with the cosmological parameters of the standard model and 
apply the algorithm to them. This would yield the probability by which such a spatial- 
correlation-function signature gives a false positive detection. However, the numerical 
effort for such an investigation is prohibitive. Thus, only ten ACDM sky maps are 
generated here using the cosmological parameters stated in Section [2j For each of the 
ten sky maps, five MCMC chains are computed and analyzed. The KQ75 mask is 
applied for all sky maps. No maxima at small values of L are found for nine out of the 
ten sky simulations where the probability ([61) is comparable or even larger than the one 
computed from the ILC map. In these cases the Markov chains are found drifting to 
ever larger values of L which means, as discussed above, a decreasing significance for 
a detection. However, one of the ten models produces such a nice peak structure as 
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Figure 11. The cross-correlation £q and the auto-correlation £4 are shown in panel 
a) for two toroidal cells: one with a correct orientation (L = 3.85) and one with a 
wrong size L = 4.00 and orientation (a = (3 — 7 = 0) using the 5yr ILC sky map. 
Panel b) displays £c and £a for the best toroidal cell for the "false positive" ACDM 
model. The KQ75 mask is applied in all cases. 

seen in Fig.HJbut around L ~ 4.17. The corresponding cross-correlation £c matches the 
auto-correlation £a even better than the best match for the torus using the real data as 
shown in Fig. [Til Furthermore, one of the remaining nine simulations produces a peak 
structure with a probability which is a factor 10 lower than that of the ILC map. Based 
on these limited observations one might conclude that ~ 10% . . . 20% of the simulations 
would produce such a false positive detection. This is a very high level and thus, the 
conclusion should rather be that, if the topology is an equilateral toroidal one, then the 
torus has probably the side length L ~ 3.85 and the orientation given in Tab.dJ 

It is intriguing that the two-point temperature correlation function C($) : = 
(5T(n) ST(n')), h • n' — cos^, also suggests a torus size around L ~ 3.86 when 
analyzed with the 3yr ILC map and applying the kpO mask |20j. Evaluating C{d) 
without applying any mask leads, however, to L ~ 4.35. But nevertheless, both the 
observed low power at large scales in the CMB sky map as well as the spatial-correlation- 
function signature are compatible with an equilateral toroidal structure. However, a firm 
conclusion cannot be drawn because of the high false positive rate. 

In [TJ the spatial-correlation-function method is applied to the Poincare 
dodecahedral space and also a favoured orientation for this topology is found. Of 
course, only one of both topologies can be realized, if one of them at all. It is currently 
not clear which topology gives a better match to the data. In a following publication 
this question will be addressed by applying the spatial-correlation-function signature to 
several distinct topologies. Furthermore, it is intended to obtain a better estimate of 
the false positive rate. 
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